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Abstract 

We study the presence of entangled states of nucleon pairs from nuclear de- 
cays and in reactions with exotic nuclei, e.g., ^^Li, or ®He. It is shown that 
the fingerprints of entangled states in these subsystems are visible in correla- 
tion measurements and can be accessed with present experimental techniques. 
This shows that not only atomic and optical systems, but also nuclear sys- 
tems serve as important tools to obtain dichotomic outcomes for tests of the 
Einstein-Podolsky-Rosen paradox. 
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I. INTRODUCTION 



In recent years, there has been much interest in the physics of unstable nuclei. Nuclei far 
from the valley of nuclear stability, with lifetime of a few miliseconds, play an important role 
in cosmology and astrophysics, as well as in the traditional applications of nuclear physics. 
Some of these nuclei possess rather unexpected properties, as for example, the two-proton 
radioactivity [1] , nuclear halos, large deformations, etc. Many theoretical models for these 
nuclei invoke subtle aspects of basic quantum mechanics. A nuclear halo, for example, is 
thought to be a simple manifestation of quantum tunneling of loosely bound states. To 
gain insight into these features, typical experiments involve the measurement of momentum 
distributions of the fragments, knock-out and stripping reactions. Coulomb excitation, etc. 

Rare nuclear isotopes can also be used to study counterintuitive aspects of quantum 
mechanics. For example, two-proton decay in s-wave states could be used for a test of 
quantum mechanics versus local realism by means of Bell's inequalities [2]. Since the final 
state of the two protons can be found in a singlet state, their wave function is spin-entangled. 
The identification of the spins of the proton in two detectors separated far away would be 
useful to test the Einstein- Podolski- Rosen (EPR) paradox [3]. In fact, these tests should 
be performed in different and complementary branches of physics to avoid the loopholes 
encountered in photon experiments. Correlation experiments with low-energy proton-proton 
scattering have already been used for studies on the EPR paradox. Two-proton radioactivity 
would also qualify for this purpose. 

Simple quantum entanglements are also visible in peripheral reactions with rare isotopes 
[4]. We will consider some examples here. To start with, we observe that in Coulomb disso- 
ciation experiments of nuclear collisions at intermediate energies 50 MeV per nucleon), 
the exchanged photon is almost real, thus the momentum and energy transferred obey the 
relation p ~ £^/c. A single nucleon cannot absorb this photon (the photon momentum is 
too small, p ~ for £^ ~ 1 MeV) . In order to preserve the energy- momentum condition the 
photon has to be absorbed by more than one nucleon. The energy can be shared by two 
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nucleons which fly apart in opposite directions. The same happens in nuclear interactions in 
peripheral collisions, when little energy is transferred to the nucleus. At forward scattering 
angles the momentum transferred is |p — p'| — p — p' cos 6' ~ p — p' ~ E /v, where p (p') 
is the initial (flnal) center of mass energy, 9 is the scattering angle, and v the projectile 
velocity. For collisions at intermediate energies, v ~ c, and the same argument as for the 
Coulomb dissociation applies. Experiments which explore this feature have been performed. 
For example. Coulomb dissociation of ^^Li projectiles (a rare isotope) into a ^Li nucleus and 
two neutrons have been done [5]. Since neither ^°Li, nor the di-neutron system are bound, 
the neutron-neutron correlation in the bound state of ^^Li is the key factor to allow the mere 
existence of this nucleus. The experiments have concentrated their efforts in understanding 
this particular aspect of the ^Li + n + n system. 

The singlet and triplet configurations of the neutron-neutron system in ^^Li have different 
contributions to the total nuclear wavefunctions and a recent experiment [6] has been able to 
access their relative weights. However, some assumptions on the reaction mechanism were 
necessary to interpret the results. It is also not clear if different spin configurations in the 
final state could change the momentum distributions of the fragments. In fact, very little is 
known about final state configurations in reactions involving unstable nuclei. 

In this article we show that one can disentangle the contributions of singlet and triplet 
spin final states of a two nucleon system in reactions with rare isotopes by measuring momen- 
tum correlations. This is a useful result as a direct measurement of the spin orientations of 
each nucleon is by far more complicated in such conditions. The application of the method is 
very general as it only relies on measured quantities, independent of models for the reaction 
mechanism. These are the widths and differential cross sections of momentum distributions, 
which have been obtained experimentally for several reactions involving unstable nuclei. 

The correlated (C) and uncorrelated (U) measurements of nucleons 1 and 2, with relative 
momentum Ap = |pi — P2I are defined by 
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where P (pi; P2) is the probabihty to measure a nucleon having momentum pi in coincidence 
with the measurement of the other nucleon having momentum p2. The integration in Qp 
is over all orientations of Ap. P (p) = / P (p; p') dp' is the probability to measure the 
momentum p for one of the nucleons, irrespective of what the momentum of the other 
nucleon is. N is the total number of particle measured, i.e., N — / P (p) dp. 
The correlation function is defined as 

One can describe each nucleon in the final state by Gaussian wavepackets: 

i/ji={27ra^) expj — j, (3) 

where the average recoil momentum of the core is Pc, and the average momentum of the 
pair at their center of mass is p„. The momentum spread of the nucleon wavefunctions 
in the final state, (T„, depends on the reaction mechanism for the specific reaction studied. 
E.g., for the breakup of ^^Li projectiles at Eiab — 300 MeV/nucleon, the neutrons acquire a 
momentum spread of the order of (T„ ~ 20 MeV/c. 

The spins of low energy nucleons imply that they can be arranged into singlet {S — 0) and 
triplet {S = 1) states. The spatial part of their wave function is symmetrized accordingly: 

^^"^^ (Pl, P2; Pc, Pn) = {V'l(Pi; Pc, Pn)V'2(P2; Pc, Pn) ± V'l(P2; Pc, Pn)V'2(Pi; Pc, Pn)} ■ 

(4) 

This wavefunction is entangled and cannot be factorized into individual particle wavefunc- 
tions. It will lead unavoidably to interferences which are visible if one observes the complete 
final state. This non-local entanglement is an example of the EPR paradox. The plus (mi- 



nus) sign refers to S = {S = 1) states, = [2(1 ± O^)] where the overlap integral i 
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given by O = / V'];(p; Pc, Pn)V'2(p; Pc, Pn)c^P = exp(-p^/8(7^). It does not depend on Pe. 
Thus, if p„ = the final wave function, ^, of the pair is 100% an 5" = state, i.e., \1' ~ 
For Pn 7^ 0, the singlet and triplet states can both contribute to "if. 



The probabilities in eq. 1 can be written as 

^(Pi;P2)= J-^^^ {|^(+)(pi,P2;Pc,Pn)f + A^|^'(-UPl,P2,;Pc,Pn)f C^PeO?p„| (5) 

where A4 is the mixing parameter, determining the relative contribution of the triplet state, 
and da/dPcdpn is the differential cross section for the process. The experiments in periph- 
eral collisions with unstable nuclei show that the momentum distributions are only slightly 
shifted from Pg = and that the shift is much smaller than the width of the recoil mo- 
mentum distribution. Thus, the correlation function is rather independent of the details of 
da/dPcdpn- However, it is important to study the dependence of the correlation function 
on p„. We thus replace Pc = in the formulas above and perform the integrals to obtain 
the corelation function i?(Ap;p„). The result can be obtained analytically. For pure singlet 
or pure triplet (-), states one gets 

^ 2 [,(:r) ±1][1± O'^] _ 

1 + 2C»4 + ^(x)C»2 ± 8/i(x)C»V2 ' 

where g — sinh(a;)/x, h — sinh(x/2)/x, and x — p„Ap/2(T^. 

The results for R^^^ (upper figure) and R^~^ (lower figure) are shown in figure 1, as a 
function of the variables Ap/an and Pn/cn- One sees that the properties of the correlation 
functions in the singlet and triplet states are completely different. When the differences of 
the average momenta of the nucleons is small Ap <^ an, the correlation function is negative 
only for the triplet state. It is -1 at Ap = for the triplet state, whereas it is close to zero 
for the singlet state. While for the former case the correlation function crosses zero at two 
points, it does not have a null point for the singlet case. 

In figure 2 we show the case for which the final state is an admixture of triplet and 
singlet states, as a function of the mixing parameter A4, which appears in eq. 5. The 
calculation also leads to closed analytical expressions. The dotted, dashed, and solid curves 
are obtained for M. — 0.1, 0.5, and 0.9, respectively. As one increases M. the fingerprint 
of triplet states in the final wavefunction becomes more and more visible. This means that 
measurements of pair correlation functions should be able to discriminate the triplet and 
singlet configurations in the final wavefunction. 



The traditional idea of diproton radioactivity is due to the pairing effect. Two protons 
form a quasiparticle (diproton) under the Coulomb barrier and this facilitates penetration. 
In a more formal description, one has a system with two valence protons in the same shell 
and coupled to J'^ — 0^. The method described above is directly applicable to determine the 
spin mixing of final states in low-energy two-proton nuclear decay for 0"*" — > 0"*" transitions. 
In this case the final spin wave function of the pair equals that of the initial wave function. 
In particular, when singlet states are identified, spin-spin coincidence experiments will gen- 
erate dichotomic outcomes for each single measurement. These can be tested using Bell's 
inequalities for spin-spin correlations along different spin-axes. Since not all the particles can 
be measured, one would have to use a softer version of the Bell's inequalities, as developed 
by Clauser and Horne [7]. In the case of peripheral nuclear reactions this is not so simple. 
But, the entanglement of the final wavefunction will also have noticeable consequences, as 
we show next. 

Let us assume for simplicity the emission of a nucleon pair in a reaction involving halo nu- 
clei. For weakly bound nucleons one can assume that the reaction suddenly detaches the pair 
from the core and that their intrinsic wavefunction is likely to be an entangled state. Their 
wavefunction can be written in the form exp(ipi.ri) exp(ip2.r2) ± exp(ip2.ri) exp(ipi.r2). 
This wavefunction leads to destructive, or constructive interferences. The cross section for 
the reaction process will be given by 

(7(P,q) = /(P,q,i?)[l±cos(gi?)]. (7) 

where q = pi — p2, P is the recoil momentum of the center-of-mass of the pair and R is 
a parameter which describes the initial spatial localization of the nucleon pair. All other 
features of the reaction mechanism are included in the function /(P, q, R). Thus, for qR -C 1 
one should be able to see a destructive interference for triplet final states and constructive 
interference for singlet final states. In fact, this was clearly seen in an experiment dedicated 
to study the existence of the nucleus [8]. The transfer reaction ^H(®He,^He)^H was 
studied by detecting two protons emitted from the decay of ^He. The energy correlation 
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of the two protons was measured, as shown in figure 3. If we apply eq. 7 to describe the 
relative energy distribution of the two protons, we get 



In figure 3 we show the result of eqs. 8 together with the experimental data of ref. 
[8]. The solid, dashed, and dotted curves correspond to singlet, triplet, and uncorrelated 
states of the proton pair, respectively. We used = 4.6 fm as the size of the source region 
where the protons originate. The agreement with the experiment is quite poor, especially 
for large relative energy of the pair. This is expected since we have neglected final state 
interactions (e.g.. Coulomb interaction). In ref. [8] this spectrum was described theoretically 
by using the Landau-Smorodinskii's cfi^cctive range approximation [9] for the protons with 
a scattering length of Upp — —7.806 fm. However, it is important to notice that the most 
important ingredient is the assumption of a singlet or of a triplet state for the protons. As 
one easily sees from figure 3 the spectrum is completely different in each situation. When 
no entanglement exists, there is no interference and the spectrum also changes dramatically. 

The examples discussed above involve nucleon pairs interacting with one or more parti- 
cles. Entanglement measures for general multiparticle systems are still under dispute [10]. 
However, if the state of the whole system is a pure state, and the full system is being regarded 
as divided into two subsystems, a convenient measure of entanglement is 



where p„ is the reduced density operator of the nucleon pair subsystem, N, and |^ > is 
the wavefunction of the total system. Under these conditions, it can be shown [11] that 
this quantity yields the same value no matter for which subsystem it is calculated, i.e., 
Pc — Pn = P (the traces in eq. 9 can be permuted). P is called the "purity" since it takes 
on its maximum value 1 if the subsystem is in a pure state. The more P deviates from 1 




singlet 



triplet 





(9) 
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the larger is the error that would occur in treating the system as an entangled state. To 
illustrate this point, let us consider the fragmentation of ^^Li projectiles, as obtained in the 
experiment [6]. 

The analysis of experiment [6] was based on Hansen's wounded wavefunction model [12]. 
This method assumes that the collision is so fast that a piece of the initial wavefunction 
is scraped-off in the collision. This allows to calculate the frozen momentum configuration 
of the initial state by a Fourier transform of this wavefunction. One therefore assumes a 
100% entanglement of the remaining constituents of the nucleus. In the experiment [6] this 
was used in the case of a knock out reaction, where a neutron was removed from ^^Li. The 
measurement of the momentum distributions of the neutron and the ^Li fragment was used 
to deduce the properties of ^'^Li. The nucleus ^^Li is usually called a Borromean nucleus, 
meaning that both the n-n, as well as the ^°Li, subsystems are not bound. We now show 
that it is exactly this property which is responsible for the successful interpretation of the 
experiment [6]. We use eq. 9 for the subsystems ^°Li (c) and n-n (N). 

We use the results of ref. [13], where it was shown that the bound-state properties of ^^Li 
can be obtained with a three-body model by using phenomenological potentials of the form 

Kn = ^ne-'''/'", K9 = S,e-''/'' + 52e-^'/^^ (10) 

where 

P = rn - r„', A = ^ (r„ + r„/) - rg, (11) 

are the Jacobian variables. The following potential parameters were used: (a) Sn — —31 
MeV, 6„ = 1.8 fm, to reproduce the low-energy scattering length and effective range of the 
n-n system, (b) The set of values Si = -7.0 MeV, Ci = 2.4 fm, ^2 = -1.0 MeV, Ca = 3.0 
fm, reproduce the size and the binding of ^^Li (~ 0.3 MeV) quite reasonably. 

In eq. 9 one should insert all bound states of the subsystems to obtain the probability 
that a measurement of the Hansen's wavefunction would give an entangled ^°Li nucleus. 
However, there are no bound states in these systems. In fact, we notice that the potentials 
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presented above barely bind the n-n and the ^Li subsystems. They yield the (negative) 
scattering lengths, a„„ — —17. 4 fm and a„9 — —10.66 fm, respectively. These are large 
values, reflecting almost bound states, with E 0. Thus, one can estimate the value of the 
purity in eq. 9 by replacing their respective wavefunction by a constant value within ^^Li. 
This yields the particularly simple expression 

( 4. < A >3 /a ) < P > Va ) "V "4/ ^' ■ '1^' 

where ^iiLi is the ^^Li wavefunction in terms of the Jacobian variables and < p >— 3.3 fm 
and < A >= 3.1 fm are the most probable values of the Jacobian coordinates. The ^^Li 
wavefunction in eq. 12 is obtained from a variational calculation, as in [13]. We get P ~ 0.8 
which is very close to unity, thus proving our assertion. 

In conclusion, we have shown that entangled states in nuclear decay and in reactions 
with exotic nuclei is another route to study important problems of basic quantum mechanics 
interest, as the Einstein- Podolski- Rosen paradox. These studies would be complementary to 
others performed in atomic physics. We have also shown that, due to the entanglement of 
final states, simple correlation functions yield compelling evidence of the spin and momentum 
states of the nuclear subsystems. The forthcoming RIA (Rare Isotope Accelerator) facihty 
will be an ideal laboratory to perform these studies [14] . 
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FIG. 1. Correlation functions for different values of the average momentum, p„. The upper 
(lower) figure is for the singlet (triplet) state. The dotted, dashed, and solid lines correspond to 
Pn/cn = 0.01, 0.1, and 0.5, respectively. 
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FIG. 2. Correlation functions for Pn/(^n = 0-5 and for different admixtures of singlet and triplet 
states. The dotted, dashed, and solid lines correspond to = 0.1, 0.5, and 0.9, respectively. M. 
is the absolute contribution of the triplet state. 
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FIG. 3. Relative energy distribution of two protons from the reaction p(®He, ppt). Calculations 
axe for triplet (T) and singlet (S) states. The dotted curve assumes no correlation in the pair 
wavefunction. 
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